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Abstract 

The solutions of the discrete Painleve equation I were constructed in terms of elliptic and 
hyperelliptic tp functions for algebraic curves of genera one and two. For the case of genus two, 
there appear higher order difference equations which naturally contain the discrete Painleve 
equation I as a special case. 

Keywords 

Discrete Painleve Equation, Elliptic ip Function, Hyperelliptic ip Function 

§1. Introduction 

In this article, we will consider the algebraic solutions of the discrete Painleve I equation 
[ORGT] 

Pn+lPn-! = A + JL (1-1) 
Pn Pn 

where z and a are some parameters. We give solutions of the equation (1-1) in terms of the 
elliptic and hyperelliptic ip functions. 

In section 2, we give an elliptic solution of this equation. Section 3 contains our main subject, 
which is based upon the recent studies on if) function [C, Ma2, MO, 01, 02]. The tp function is 
defended over an algebraic curve itself embedded in its Jacobian rather than over the Jacobian 
variety. Although they need slightly corrections, Cantor essentially gave a determinant expres- 
sion of ^-function and a recursion equation on the ip- functions of genus two [C]. On the other 
hand, Onishi gave another determinant expression of ip- function [02]. Recently both expressions 
are connected by us [MO]. In §3, we show that the recursion relation with a correction become 
(1-1) and a natural third order difference equation for a certain point in the related algebraic 
curve, 

4+24-1 = , " 5 , "4 [-7- 1" ~f - ) 1 I 1 " 2 ) 

where a's are some parameters. (1-2) should be regarded as an extension of (1-1). By using 
our identification [MO] and the expressions of Onishi, we give their solutions in terms of the 
hyperelliptic a functions [Bl-3]. We also obtain the sixth and forth order ordinary difference 
equations related to genus two curves. 
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This study started after listing to the lecture of the recent progress of third order difference 
equation by Yahagi, Kimura, Tusjimoto and Hirota [YKTH, HYK]. They have been studying 
the third order difference equations. Their motivation is to construct a list of the integrable 
ordinary third order differential equations as Painleve school gave one for second order one in 
beginning of the last century. Recent studies of difference equations show that integrablity, at 
least, for the case of the ordinary second differential equations is determined by properties of 
their difference equations versions. Hence in order to construct the list, they started to classify 
the third order difference equations by aided of numerical computations. They have a list which 
contains eleven types of the third order difference equations. They stated that they found an 
equation whose solution can not be expressed by any elliptic functions. As 1 believe that some 
class of the difference equations must be defined over an algebraic curve itself embedded in the 
Jacobi variety whereas the continuous soliton equation is related to theory of Jacobi variety 
rather than curve itself. Thus 1 began to study this problem along the line of arguments of 
Onishi [01,2] and Cantor [C]. 



2. Elliptic V-Function 



In this section, we consider an elliptic curve, 



y 2 =4x 3 -g 2 x-g 3 , 



(2-1) 



where g's are complex numbers. 

The elliptic V'-function [W] is defined as 



V'nH = 



a(nu) 
a(u) n2 



(2-2) 



where cr(u) is the Weierstrass sigma function. This function has Brioschi-Kiepert formula [W-W, 
p. 460 and refecences in 02], 



^(u) = (-l)™(™- 1 )/ 2 (l!2!---n!) 



2 a (. nu ) 
a(u) n2 



p'(u) p"{u) 
p"(u) p"'(u) 



The ijj function also obeys the recursion relation, 



p^(u) 
p (2 — 3 >(u) 



(2-3) 



(2-4) 



This relation is proved by the additive formula of Weierstrass p function [T] . We note that this 
recursion relation differs from the identities of Hankel determinant (2-3). For m = 2 case, we 
have a bilinear difference equation, 

1p n+ 2tpn-2 - 1p%1pn-llpn+l + ^Z^l^n^n = 0. (2-5) 

By introducing the quantity (3 n = V'n+iV'n-i/V'n' (2-5) turns out to be the difference equation, 



Pn (3; 



2 ' 



(2-6) 



This equation is identified with the discrete Painleve I (1-1) [ORGT]. We have a special solution 
of (1-1) as 

_ cr({n + l)u)a((n - l)u) 

Pn — , no / \o ■ \ Z ' ) 



}3. Hyperelliptic ^-function of Genus Two 



In this section, we will deal with the ^-function over a hyperelliptic curve C of genus two 
defined by an affine equations, 



A + Xix + X 2 x 2 + A 3 x 3 + A 4 x 4 + A 5 x 5 , 



(3-1) 



where A5 = 1 and Aj's are complex numbers. We denote its corresponding Jacobi variety by J 
and image of inclusion of the curve into J by l(C). 

The hyperelliptic ip function of genus two is given by [01, 02], 



ipn(u) 



a(nu) 
<y 2 {u) n2 



Here u := (u\,u 2 ) is a coordinate restricted to curve itself i(C), 

^ x ' y) dx f {x ' v) xdx 



ui 



j 00 y J 00 



v 



in the Jacobi variety J := {(ni,U2)}, 



{x ' y) xdx f {x2 ' y2) xdx 



Ui 



r (x ' y > dx + r (X2 ' V2 > dx _ n x ' y > xdx | r 
j 00 y j 00 y j 00 y Jo 



V 



(3-2) 



(3-3) 



(3-4) 



for (x,y), (^2,2/2) G C. Further a is Baker's sigma function [Bal, Ba2, Ba3] and a 2 is its 
derivative with respect to u 2 . In this article, we will follow the arguments in [01, 02, MO]. 
We note that u\ is a function of u 2 [MO]. This ip n (x,y) is a polynomial whose zero (xo,yo) 
is necessary and sufficient condition for the element r ■ (x ,2/o) of J to lie in i(C) again. We 
emphasize that the function ip n is defined over the curve itself t(C) rather than J . This definition 
is key of these studies and was, first, given by Grant [G]. 

Cantor showed that the ^-function can be expressed in terms of Hankel determinant [C]. 
However his expression slightly needs a correction on its factor on y and in [MO, Ma2], we 
will give a correction along the line of argument of Onishi [02, MO]. (For example, if) 2 of 
genus two must be proportional to y from the argument in [01] but in [C] V2 is constant.) 
In our new derivation of ^'s [MO], the Toelpliz determinant is more natural than the Hankel 
determinant (see Appendix). Thus in this article, instead of Hankel determinant, we use a 
Toelpliz determinant, 



g[m+n— 1] g[m+n-2] 
g[m+n+l] g[m+n— 1] 



g[m+2n-3] g[m + 2n-4] 
[m+2n-l] g[m+2n-3] 



j[m+l] 
Jm+2] 



Jm+1] 



g[m+n — 1] g[m+n — 2] 
g[m+n] g[m+n— 1] 



(3-5) 



and T : 
s and 



i-n [ g(s), -j- ) = where m and n are positive integers, g(s) is a function of an argument 



1 d v 



-9(s). 



n\ ds n * 

As we showed in [MO], we have an expression of if) n as (see Appendix), 



(3-6) 



y 



^- 1)/2 -T^ 3)/2 (y,£) fornodd 



{n - 1)/2 -T%l 2)/2 (y,£) forneven 



(3-7) 



Noting y 2 = f(x), y 2n 1 d n y/dx n is a polynomial of x and coprime to f(x) in general. Hence 

y n(2m+2n-3) T (m) ^ Qr 



y2m+2n — 3y[m+n — l] y2m+2n — 5y [m+n — 2] 



( 2m+2n— ly [m+n] 



y2m+2n-3y[m+n — 1] 



y2m+ly[m + l] 
y2m+3y[m+2] 



y 



y2m—ly[m] 
2m+l [m+1] 



y2m-\-4n — 7y [m+2n — 3] ^2m+4n — 9^ [m+2n — 4] # , , y2m-\-2n — 3y[m-\-n — l] y2m-\-2n — 5y[m-\-n — 2] 



y2m+4n-3y[m+2n-2] y2m+4n-7 y[m+2n-3] 



y2m+2n+l y[m+n+l\ y2m+2n-3 y[m+n- 1] 



is an element of C[x] and coprime to y 2 . Hence tp n (u) can be expressed by 



8y 3 a n (x) for n = odd 
2ya n (x) for n = even 



(3-8) 
(3-9) 



where a n (x) is a polynomial of x and coprime of y. 

Cantor, essentially, showed that this ip n obeys a recursion relation, which is an extension of 
(2-4) [C, Ma2], 



0>20m0nV>n+m0 m — n 



■ipm-2lpn 0m-10n+l 0m0n+2 
0m-10n-l Ipm^n 0m+10n+l 

"0m0n-2 m +10„-l 0m+20n 



(3-10) 



(3-10) slightly differs from original one in [C] because original one needs a correction [Ma2, MO]. 
We emphasize that this relation (3-10) characterizes the ip function over a genus two curve and 
cannot be reduced to some identities of elliptic functions in general because (3-10) holds due to 
the addition relations of genus two [Ma2] . 

Noting 0o = 0, 0i = 0, tp2 = 2y and 03 = 8y 3 , we have n = 3 case of (3-11), 



0302 - 0m+3'0m-3'0m ~ 3 $ -0m+20m0m-2 ~ -0502 VVn+lVVnVVn-l 
+ 04^302 (Y>m-20m+l + 0m+20m-l) = 0, 



(3-11) 



and n = 4 case, 



-04-0 2 0m+40m-40m - 3 = Vm+20m0>m-2 + 06 "04 "02-0™ ~ (0502 + 0603 )0m+10m0m-l 

+ 050403 (0m-20™+l + 0m+20i-l) = 0. (3-12) 
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(3-11) and (3-12) are reduced to [C], 

^4^2^71+4^71-4 ~ ^5^3^2^71+3^71-3 + (^5^3 ~ ^4^)^+2^-2 

+ folplfolpn+llpn-l ~ IpeM^nlpn = 0. (3-13) 

This bilinear difference equation (3-13) is a genus two analog of the bilinear equation (2-5) and 
also characterizes the ip function over a genus two curve. We emphasize that from point of 
view of bilinear difference equations, (3-13) is a very natural extension of (2-5) and can not be 
reduced to a three-term relation in general. 

From here, we will consider modifications of (3-10)-(3-13) to find difference equations which 
have form, 

fl(x n +m-l, ,X n - m ' + l) (-IT A\ 

■^n+m^n — m' » , \i V° ± ^ t / 

J2\X n +m— lj ' i %n — m' + lj 

where /'s are polynomials of x n+m -i, • • • , x„_ m '_|_i and m, m' = 1, 2, 3, • • • , by deforming (3-13) 
as a generalization of (2-6). 

By setting b n = i/) n+ 2ipn-2/ipni (3-13) becomes a sixth order ordinary difference equation, 

MlK+sK-S ~ Jp p ^3 -gr~ (^3^6n-2«£_l&n&n+l ft n+2 

- (V'sVl - iplip2)b n -ib 2 n b n+ i + ipefpl&K ~ ^6^4^/1) = 0, 

(3-15) 

for a general point in f-(C). (3-15) can be regarded as a generalization of (2-6). 

Further we will consider other difference equations by dealing with special points. First we 
will deal with points (x, 0), y = 0. For the point, (3-13) becomes 

a n+4 a n _ 4 - (al)a n+2 a n - 2 + a 6 a n a n = 0. (3-16) 

When we define c n = a n+ 2a n -2 / 'a^, (3-16) turns out to be the discrete Painleve I (1-1), 

c n+2 c n _ 2 = j". (3-17) 

In other words, we have a hyperelliptic function solution of the discrete Painleve I [ORGT] (1-1) 
again, 

_ g((B + 2y(n-2) M ) 

Cr(ntt) 2 (T2(ti) 8 

where u is 0)). 

Further using (3-16), we have that for a point u G i((x,0)), 

d . = ip m +i(u)ip m -2(u) = a((m + l)u)a{(m - 2)u) 
m ' ij} m {u)ij) m -i{u) a{mu)o{{m - l)n)cr 2 (n) 4 

obeys the third order difference equation, 

d m+2 d m _i - - " 5 - + a 4 ( — + ^- ) = 0. (3-20) 

«m+l«m \"m+l «m / 
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We note that although n vanishes for a point i((x, 0)), d m has non-trivial values because it can 
be expressed by 

:= g^^ , (3-21) 
a m (u)a m -i(u) 

In other words, the third order difference equation (3-21) is well-defined and its solution is 
associated with an algebraic curve. (3-21) can be also regarded as a generalization of (2-6). 

Next we will deal with a point satisfying 4 = 0. Then we apparently obtain a forth order 
ordinary difference equation, 

050302 6 n+2&n-2 + 72 TTTo (0503 & n-l&n fe n+l ~ 0603 02&n) = 0, (3-22) 

°n-l n°n+l V 7 

at the point. However if 4 vanishes at a point P and P is a forth cyclic point, it means that 
02 = 03 = 04 = 0, 03 = 04 = 05 = or 4 = 05 = 06 = from the definition of function 
of genus two [C,01]. Thus after fixing the curve (3-1), we need more precise consideration on 
(3-22) whether it has non-trivial meaning or not. 



54. Discussion 



As mentioned in the introduction, the authors in [HYK] and [YKTH] have a list of ordinary 
third order difference equations. Our third order difference equation (3-20) also appears in it. 
Yahagi shows [Y] that (3-20) is a special case of (86) in [HYK], 

a + a x x n + aix n+ i + a 3 x n x n+1 

x m+2 x m -i = — — — . (4-1) 

fl3 + bix n + OlX n+ i + b 3 x n x n+1 

where a's and 6's are parameters. (4-1) becomes (3-20) by letting a 3 = bi = 0, a\ = —a^ 
and a = q 5 . Though I started this study to find a hyperelliptic function solutions of the 
third order difference equations whose solutions can not be expressed by elliptic functions, our 
attempt failed. As n at i((x, 0)) becomes the discrete Painleve equation I, which has the elliptic 
function solution, (3-20) is reduced to the ordinary second order difference equation [YKTH, Y]. 

However (3-20) is an equation at the special points satisfying y = 0. We have the bilinear 
equation (3-13) as an extension of (2-5) and (3-10) as that of (2-4) and (3-10), (3-13), (3-15) 
(and (3-22)) cannot be expressed by some identities of elliptic functions in general. Even though 
they cannot be reduced to the third order difference equation expect (3-20), it is expected that 
they might be connected with higher integrable difference equations. 

Further it is also expected that the recursion relations of 0-functions in more general algebraic 
curves, if exist, might contains interesting several difference equations. In fact, even though it 
is not correct by primitive considerations, a naive extension of the right hand side of (2-4) and 
(3-10) might be given by 



0m-30n m -20n+l m -10„+2 0m0n+3 

0m-20n-l 0m-10n 0m0n+l m+ 10 n +2 

m _10 n _ 2 m _l0 n _l m +10n m+2 0n+l 

m n _ 3 m 0n-2 m+2 0n-l 0m+30n 
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(4-2) 



which looks to contains various third order difference operator. Thus I believe that it is important 
to study the recursion relations of ip functions over more general algebraic curves [C, Ma2, MO, 
03, 04]. 

Next we will comment on a relation of our theory to so-called Sato theory [SN, SS]. As shown 
in [MO], in the derivation of (3-7) we encounter the Wronskian. We will give a rough sketch 
of the derivation of the proof in the Appendix. On the other hand in the Sato theory, the 
Wronskian also plays the important roles [SN, SS]. In the theory, we encounter a differential 
equation, 

dW m 



8t r 



B n -W m -W m -d n , (4-3) 



where W m = d m + wid m ~ l + • • • + w m , d := d/dt\, w a are functions of the Sato coordinate 
t = (t i , *2 j ■ ■ ■ ) s an d B n is a certain differential operator [SN]. For the independent solutions 
{fi} of W m {ti = 0, i > l)f(ti) = 0, we can define a Wronskian, 

7»(t) = det^Dx^D^/^ti)). (4-4) 

Roughly speaking this Wronskian becomes r function, which plays central roles in modern soliton 
theory, r is defined over a Jacobi variety of corresponding algebraic curves [Mul] . In the theory, 
algebraic curve is an auxiliary object after we obtain the Jacobian. As the Jacobian is realized 
as a quotient complex vector space divided by a discrete lattice, the r function is also related to 
a discrete equation of a formula of Jacobi variety J , which is known as Fay trisecant formula 
[F]- 

On the other hand, our theory is of functions over a curve i(C) C J based on theory of 
hyperelliptic sigma function [B1-B3]. As I pointed out in [Mai], in Sato theory we consider 
the behavior of the differential form of the second kind around the infinite point, whereas in 
Baker's theory, the differential of the first kind, which is homomorphic all over the curve, plays 
central roles. The r function can not be explicitly expressed in general except soliton solutions 
or elliptic function solutions, but in Baker's theory all quantities concretely can be expressed 
any points on a curve [Mai, BEL]. The Baker's sigma function theory and r function theory are 
resemble but slightly different. 

Further in general, functions defined over Jacobi variety, i.e., r and a functions, do not behave 
well over a curve itself; they vanish or diverge at t(C). Using the properties of a functions, we 
can tune the functions and define ip functions over there as in (3-10) [01, 02, Ma2]. In our 
theory, important quantities are defined over a curve itself. For example, the Wronskian (A-10) 
appearing in the derivation of (3-5) in [MO] is a function over a curve except infinite point. 
Thus r and ip functions are also different. Similarly I think that (3-11) can not be obtained by 
a ordinary modification of Fay trisecant formula [F]. In other words, I believe that our theory 
has an aspect which one cannot reach by using Sato theory. 
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§Appendix Determinat Expressions of ^-functions 



In this appendix, we will roughly review of the dterminant expression of ^-functions for genus 
two along the line of the arguments in [MO] in order to show that the ^-function is related to 
Wronskian, which is similar to the situation of Sato theory [SN, SS]. 

As ^-function is a polynomial in x and y over the curve C (3-1), the polynomial belongs to 
an algebra R = C[x,y]/(y 2 — f(x)). For the zeros (xj,yj) (i = 1, • ■ ■ , n 2 — 1) of ip n (x,y), there 
exists a point P in the curve C, n • (xi,yi) belongs to l(c), n-(xi,yi) — n-oo = P — oo in the sense 
of Picard group Pic°(C). In other words, for a local parameter u 2 of t(C) C J, u = (111(112), 112), 
there exist functions f,heR over the curve satisfying, 

f(v 2 ) = v2h(v2), (A-l) 

or 

^/k=o = 0, 0<m<n, (A-2) 

for V2 = U2 — u 2 ((xi,yi)), where d U2 := (d/du 2 )\u 2 - I have employed the notations d U2 as a 
partial differential following the notations of theory of a complex manifold. However it should 
be noted that u 2 and v 2 are local parameter of the curve C, and u\ is a function of n 2 . Thus 
we can replace d U2 with the ordinary differential operator (d/du 2 ) as we regard R as a set of 
functions only of u 2 (and t> 2 ). Further we notice that as J can be identified with Picard group 
Pic°(C), the addition in this appendix is defined in J; generally we cannot define addition in 
i(C) and C itself. 

The elements in R has a natural order with respect to the order of divergence at the infinity 
point, which is a natural extensions of the degrees of polynomial. We regard R as C-vector 
space generated by an ordered set, R =< 1, x, x 2 , y, x 3 , xy, x 4 , x 2 y, ■ ■ ■ >. By introducing new 
notations, 

{x s , for m < 2 

x p , for p = (s + 2)/2, m > 2, m is even , (A-3) 
x 9 y, for q = (s — 3)/2, m > 2, m is odd 

R =< l,ifi,(f2, - ■ ■ ■ >■ Let us consider its one- forms, 

=< d(p!,d(p2, ■ ■ ■ ,d(p n -t > . (A-4) 

Each c?99 m is expanded by, 

9 2 d 3 
d<fm = q^P™ ' du2 + ^ ^™ ' (du 2 )' s>2 + 7^3 • (du 2 )® 3 H • (A-5) 

We introduce its truncated expansion, 
B r) 2 

d^ ■■= • du 2 + • (du 2 r 2 + ■■■ + TT^T^m ■ (du 2 r {n - 1 l (A-6) 

OU2 ou 2 ou 2 

and a (n — l)-dimensional C-linear space, 

1/W =< dn 2 , (du 2 )® 2 , • • • , (dn,)®^" 1 ) > . (A-7) 



Let us consider the C-linear map at each point u, 

. V (n) _^ dR (n)_ ( A _g) 

At the zero of ip n , {xi,Vi), (A-2) implies that there exists a polynomial (fi^ := ^" =s a ra (f m 
(a m G C), 

q^<P H \u 2 =u 2 ( x ^)=0, 0<m<n. (A-9) 

It means that d0( n ) := X^"=m a m<kPm vanishes at «2 = U2(xi,yi)- In other words, the rank of 
the map is smaller than n — 1 and the Jacobian of 

J {n \u) = det (n _ 1)x(n _ 1) (A-10) 

vanishes. This (u) can be regarded as the Wronskian and as a coefficient of a natural form 
(du 2 )® k, s, W^ n \u) := J(")(«) • (dns)®"^" 1 )/ 2 . 

Thought J(") (u) is proportional to ^> n , J ( - n - ) ( n ) i s n °t polynomial in x,y in general. Thus we 

d d 

will modify it using a formula — — = — — in J and a relation over t(C), 

OU2 x ox 

( d j r A /«\«(«-i)/2 / 
det (n _D x(n _D I — ^ ( l) I = (-) det^.Dx^.D (^¥> W J • (A-ll) 

We note that as R is a set of functions of x (not of x), we can also replace <9 J /cte J with an 
ordinary differential operator d? jdxK It turns out that x ( ~ 9 ~ 1 ~ >n( - n ~ 1 ^' 2 (u) is a polynomial of 
(x,y). Accordingly we define 

ij) n {u) =7 /(«-i)/2jW( M ), (A-12) 

where 7 is a constant factor, which will be determined later. 

Here we permute the generator of R in terms of a permutation group 6, 

dR {n) = &dR (n) =< dx, d(x 2 ), ■■■ , d(x p ),dy, d(xy), • • • , d(x q y) > 

=:< dipi,dip 2 , ■ ■ ■ ,d<p n -i > . (A-13) 

where p G N and q G Z> ; 

_ j n — 1, for n < 3, 

P ~ I > for n > 3 > 

not defined, for n < 3, , . , . N 

[^], for n>3. (A-14) 



Table 1: 5 = 2 

r |123456789 

p 1012233445 
g j - - - 1 1 2 2 

p-g+1 j - - - 3 4 3 4 3 4 
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The quantity det( n _i) X ( n _i)(5|^w) is computed by 

rP -i 



J(n) 



1 2x 
2 



px h 



U 



p(p — l)x p 2 y 



(xy) 
(xy) 



p \ y(p) ( X y)(p) 

y (p+l) ( X y)(P+l) 



,{n-l) ( xy )(n-l) 



(xiy) 
{x'iy) 

( x q y )(p+l) 

(x1y)(n-l) 



3! 



pi 



,(p+l) 



,{n-2) 
Jn-l) 



„_ 2 P l2/ ("- 3 ) 



-iP q -iy<*-«> p+1 P q y( p -i+V 



q—ly 

„P , ? /(p) 
n—2 1 q—ly 



„P ? /(p) 
n — 2 1 qy 

„ P 7/CP+ 1 ) 
n — 2 1 qy 



1!2! • • -p!(p + l)!(p + 2)! • • • (n - 1)! • T^- q+1) (y, A), 



(A-15) 



where n P m := n!/(n — m)!. Let 7 = det(S)/l!2! • • • (n — 1)! in (A-12) and then we have 

d 



lPn{u) =yMn-l)/2 T (n-l ){y a y 



9+1 

This is reduce to (3-7) and ?/> n can be regarded as a Wronskian. 



(A-16) 
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